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Abstract 

In this note, we examine the problem of identifying the interaction geometry among a known number 
of agents, adopting a consensus-type algorithm for their coordination. The proposed identification process 
is facilitated by introducing "ports" for stimulating a subset of network vertices via an appropriately 
defined interface and observing the network's response at another set of vertices. It is first noted that 
under the assumption of controllability and observability of corresponding steered-and-observed network, 
the proposed procedure identifies a number of important features of the network using the spectrum 
of the graph Laplacian. We then proceed to use degree-based graph reconstruction methods to propose 
a sieve method for further characterization of the underlying network. An example demonstrates the 
application of the proposed method. 

Keywords: Inverse problems, coordination algorithms, system identification, graph reconstruction. 

I. Introduction 

Physical sciences are often concerned with inferring models and physical parameters from 
data. Given a model for a physical phenomena, computing the data values is often referred to as 
the forward problem. On the other hand, in inverse problems, the objective is the construction, 
validation, invalidation, or reconstruction of the model from a set of measurements associated 
with the system. Inverse problems arise in fields such as astronomy, geophysics, medical imaging, 
remote sensing, ocean acoustic tomography, and non-destructive testing [1], [2], [3]. Closer to the 
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present work are the inverse problems associated with electrical networks [4], and the celebrated 
"Can one hear the shape of a drum?" which aims to characterize a manifold via its spectra [5], 
or more recently, "Can one hear the shape of a graph?" [6]. In fact, in this paper, we address the 
inverse problem related to consensus-type coordination algorithms. Consensus-type algorithms 
have recently been employed for analysis and synthesis of a host of distributed protocols and 
control strategies in multi-agent systems, including, flocking, formation control, rendezvous, and 
distributed estimation [7]. 

One of the key aspects of this class of protocols is the strong dependency between the 
interaction and information-exchange geometry among the multiple agents, on one hand, and the 
dynamic properties that these systems exhibit, on the other. Motivated by this dependency, in our 
work, we consider the scenario where the interaction network is inside a "black box," and that 
only certain "boundary" nodes in the network can be influenced and subsequently observed. The 
"input" boundary nodes are then used to stimulate the network, whose response is subsequently 
observed at the "output" boundary nodes. Using this setup, in our complementary work [8], we 
have presented a node knockout procedure that aims to find the generating function of the graph 
Laplacian from the observed input-output data. Our focus in the present work, in the meantime, is 
to reduce the search space for the identification of the network topology by blending ideas from 
system identification, integer partitioning, and degree-based graph reconstruction. The implicit 
contribution of our analysis is its ramifications for exact identification from boundary nodes 
for networks that have an embedded consensus-type algorithms for their operation, including 
formation flying, distributed estimation, and mobile robotics. 

Our notation and terminology are standard.' We denote hy Q = (y,£) the undirected simple 
graph with vertex set V and edge set £, comprised of two-element subsets of V; we use "nodes" or 
"agents" interchangeably with "vertices." Two vertices u,v e G are called adjacent if {u, v} e S. 
For vertex i, deg i denotes the number of its adjacent vertices or neighbors. The Laplacian matrix 
for the graph Q is denoted by L(^). Laplacian matrices are positive semi-definite whose spectrum 
will be ordered as = Ai(L(^?)) < X2{L{g)) < . . . < An(L(^?)). We use (f)g{s) to denote the 

'The main focus of this work is on undirected graphs. However, the extension of some results to the directed case will be 
examined in subsequent works. 
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characteristic polynomial of the graph Laplacian. The cardinality of the set l-L will be denoted 
by \H\; 0{f{n)) on the other hand denotes a function of n that is bounded by some constant 
multiple of f{n) for large values of n. 

II. Network Identification 

Consider the consensus protocol adopted by n-nodes, where Xi is the state of the i-th node, 
e.g., its position, speed, heading, voltage, etc., evolves according to the sum of the differences 
between the i-th node's state and its neighbors. Next, let a group of agents X c V with cardinality 
|X| = rj, "excite" the underlying coordination protocol by injecting signals to the network, with 

another set of agents O G V, of cardinality \0\ = ro, measuring the corresponding network 
response. Hence, the original consensus protocol from node i's perspective assumes the form 

ii{t)= ^ {x,{t)-x,{t)) + B,u,{t), (1) 

where Bi — /3i if i e X, and zero otherwise. Without loss of generality, we can always assume 
that /5i = 1 and modify the control signal Ui{t) as l3iUi{t) if necessary. Adding the observation 
ports to this "steered" consensus, and having yj{t) = Xj{t) when j e O, we arrive at the compact 
form of an input-output linear time-invariant system, 

x{t) = A{g)x{t) + Bu{t), y{t) = Cx{t), (2) 

where A{g) = -L{g) e R"''", B e R"^''^, and C e R'^o^". 

Even though in general sets X and O can be distinct and contain more than one element, 
for the convenience of our presentation, we will assume that they are identical- and at times, 
assume that the resulting input-output system is in fact SISO. The extension of the presented 
results to the case when X and O are distinct will be discussed after introducing the basic setup 
and approach. 

We now pose the inverse problem of graph-based coordination algorithms, namely, the feasi- 
bility of identifying the spectral and structural properties of the underlying network Q via the 
data facilitated by the input-output ports X and O. In order to implement this program, however, 
we need to assume that: (1) the identification procedure has knowledge of the number of agents 
in the network, and (2) the input/output sets X and O have been chosen such that the system 
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described in (2) is controllable and observable. Although the first assumption is reasonable in 
general, the second one requires more justification which we now provide. In the trivial case when 
X = V and B is equal to the identity matrix, the input-output consensus (2) is controllable and 
by duality, observable. However, more generally, the controllability /observability of the network 
from a subset of its boundary nodes, is less trivial, and more to the point, not guaranteed for 
general graphs [7]. In the meantime, since we will need controllability and observability of the 
network for its identifiability, we will rely on a topical conjecture in the algebraic graph theory 
community to the effect that for large values of n, the ratio of graphs with n nodes that are 
not controllable from any single node to the total number of graphs on n nodes approaches 
zero as n — )> oo [9]; this phenomena is depicted in Fig. 1. In the present paper, we take the 
controllability and the observability of the underlying graph from the input and output nodes as 
our working assumption. In the meantime it is always convenient to know when the network is 
uncontrollable from a given node. 

Lemma II.l. Let G{s) — C{sl — A)~^B as the input-output realization of (2). The uncon- 
trollable/unobservable eigenvalues of (2) will not appear in the corresponding entry of G{s). 
Specifically, G{s) will be order n — i polynomial for the SISO case with n agents and i 
uncontrollable/unobservable eigenvalues. 

Proof: Since the underlying graph is undirected, matrix A{Q) is symmetric and there exists 
a unitary matrix U and a real nonnegative diagonal matrix A = diag(Ai, . . . , A„) such that 

A. {<3) = UAU^. The columns of U are an orthonormal set of eigenvectors for A{Q). The 
corresponding diagonal entries of A are the eigenvalues of A(Q) [10]. Therefore, 

G{s) = C{sl - A{g))-^B = C{sl - UAU^y^B = CU{sI - A)-^U^B (3) 

From PBH test, if the system (2) is not controllable, there is an eigenvector that is orthogonal to 

B. Therefore for an arbitrary uncontrollable eigenvalue Aj, the i-th row of U'^ is orthogonal to B 
and Aj will not appear in {si — A)~^U'^B. An analogous argument works for the unobservable 
case. ■ 
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A. System Identification 

We now consider various standard system identification procedures in the context of identi- 
fying the spectra of the underlying graph Laplacian, and subsequently, gaining insights into the 
interconnection structure that underscores the agents' coordinated behavior. 

System identification methods are implemented via sampling of the system (2) at discrete time 
instances,^ 5, 25, ... , kS, . . ., with S > 0, assuming the form 

z{k + l)^Adz{k) + Bdv{k), w{k)^Cdz{k), (4) 

where z{k) = x{k6), v{k) = u{k5), w{k) = y{k6). Ad = e^^, Ba = e^'dtj B, and Cd = C? 
In fact, the system identification process leads to a realization of the model 

z{k + 1) = Adz{k) + Bdu{k), w{k) = Cdz{k), (5) 

where [Ad, Bd, Cd) is the realization of [Ad, Bd, Cd) in (4). The estimated system (5), on the 
other hand, is equivalent to the continuous-time system 

x{t) = Ax{t) + Bu{k), y{t) = Cx{t), (6) 

with Ad — e^-^, Bd—(^ e^^dtj B, and Cd — C] in this case. A — [1/5) logj^ Ad where \og^ 
denotes the matrix logarithm. Since the system (5) is a realization of the system (4), it follows 
that the estimated triplet {A, B, C) is a realization of {A, B, C) in (2). As a result, there exists 
a similarity transformation induced by the matrix T, such that A = TAT~^, B = TB, and 
C = CT~^. In fact, in the controllable/observable case, eigenvalues of Ad are precisely matched 
with the eigenvalues of Ad. Obtaining a zero as eigenvalue of Ad, which is equivalent of obtaining 
— oo as the eigenvalue of A, is a sign of uncontrollable and/or unobservable mode in {2)^ For 
example in the identification procedure called Iterative Prediction-Error Minimization Method, 
the model (4) for every input vi and output Wj can be represented as A{q)wj{k) = 'B{q)vi{k), 
where A{q) = 1 + aiq~^ H h a„g~" and B(g) = biq~^ H h hrj^q~'^^. The unknown model 

^The system identification methods work based on data sampling from tlie system. Since we aimed to identify the interaction 
geometry of the network, we originally considered a continuous system. Therefore, we need to discretize the system (2). 
^The notation for a square matrix A refers to its matrix exponential. 

''This follows from Lemma n.l since — oo will appear as zero in the corresponding entries of P{s). 
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parameters 9 — [oi, . . . , a„, 61, . . . , fe^.^.] can then be estimated by comparing the actual output 
Wj{k) and predicted output Wji{k\k — 1) using the mean-square minimization. In this case, the 
output predictor is constructed as Wji{k\k — \) = [—Wj{k — n), . . . , —Wj{l),Vi{k — rx), ■ ■ ■ ,Vi{l)]. 
In yet another candidate system identification procedure, namely the Subspace Identification 
Method, the system (4) is approximated by another system in the form (5), using the state 
trajectory of the dynamic system that has been determined from input-output observations. The 
Hankel matrix, which can be constructed from the gathered input-output data, plays an important 
role in this method. By constructing the Hankel matrix, the discrete time system matrices A^,, 
Bd, and Cd can then be determined. Subsequently, the continuous-time estimated matrices A, B, 
and C can be identified; see [11] for an extensive treatment of system identification methods. 

In summary, an identification procedure such as the above two methods, implemented on a 
controllable and observable steered-and-observed coordination protocol (2), leads to a system 
realization whose state matrix is similar to the underlying graph Laplacian and in particular 
sharing the same spectra and characteristic polynomial. However, a distinct and fundamental 
issue in our setup is that having found a matrix that is "similar" to the Laplacian of a network is 
far from having exact knowledge of the network structure itself [8]. This observation motivates 
the following question: to what extend does the knowledge of the spectra of the graph, combined 
with the knowledge of the input-output matrices, reduce the search space for the underlying 
interaction geometry? In this note, we explore this question using techniques based on integer 
partitioning and degree-based graph reconstruction. 

III. Graph Characterization 

We first review qualitative characterization of the underlying interconnection topology via its 
identified characteristic polynomial. We then explore the possibility of reducing the search space 
for the underlying network via the proposed sieve method. 

A. Graph Characterization via the Characteristic Polynomial 

Recall that via a system identification method, the characteristic equation of the system (2) 
can be found as 

^g{s) = dei{sl - A{g))^ + ais^-^ + ... + an-is + an. (7) 
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Although the spectra of the graph Laplacian in general is insufficient to form an explicit 
characterization of the underlying network, it leads to a number of useful structural information 
about its geometry; we list a few: (1) the value (1/?^) nr=2 ^*(^) number of spanning 

trees in Q, (2) one has \S\ = 1/2 X^"^^ Aj, where \£\ is the number of edges in the graph, (3) 
if a„_i = n, the underlying interconnection is a tree. For a tree, the coefficient a„_2 is also 
called the graph Wiener index (the sum of all distances between distinct vertices of Q) [12], (4) 
if the associated graph is a tree, is the number of A; -matching in the subdivision of Q (where 
each edge of Q is replaced by a path of length 2), (5) if the eigenvalues of L(Q) are distinct, 
with 0, A2, . . . , \r where A,. > . . . > A2 > 0, define V'e(^) = (x — A2) . . . (,t — X,). It is then 
well-known that ijjgo = (—lYtjjg(n — x), where the graph Q"^ is the complement of the graph Q. 
The Hoffman number of the graph, can also be found as X2X3 ... Xr/n whose properties 
and applications have been studied in [13], (6) let T be a tree with n > 2 vertices. If T has 
only one positive Laplacian eigenvalue with multiplicity one, then T is the star -fri,„_i [13], (7) 
let ^ be a connected graph with exactly three distinct Laplacian eigenvalues. Then the algebraic 
connectivity (the second smallest Laplacian eigenvalue) of is equal to one if and only if Q is 
a star of Ki n-i with n > 3, (8) if is a connected graph with integer Laplacian spectra, then 
d{Q) < 2k(Q), where d{Q) is the diameter of the graph and k{Q) = (1/n) 111=2 

Although the spectra of the Laplacian provides important insights, as noted above, into the 
structural properties of the network, we now proceed to explore the possibility of complete 
identification of the underlying network using its graph spectra complemented with a sieve 
method. 

B. Graph Sieve 

In this section, we provide an overview of the graph sieve procedure- that in conjunction 
with the identified Laplacian spectra- leads to a more confined search for the network in the 
black box. The essential ideas involve the judicious use of integer partitioning algorithms and 
degree-based graph reconstruction. 

Recall that with the standing assumption of C = in (2), for the identified system matrices 
(C, A, B), after appropriate relabeling, the product CAB = CAB leads to the first r x r block 
partition of the matrix L{Q). Notice that if B ^ C^, we still obtain entries of the matrix 
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L{Q) which may not contain the diagonal entries. The product CAB gives us the degree of the 
nodes that are in common in both input and output sets and information on the minimum degree 
of other nodes. Since the eigenvalues of the identified matrix A are identical to those of L(Q), 
as the result of the identification process, we have access to the sum of the degrees of all nodes 
in the network, as well as the degrees of a subset of r-boundary nodes. Let us define R as the 
set of nodes in r-boundary nodes which appears in both I, and O with \R\ —f. Moreover, let 



If _B = C^, then = iriice{C AB), and the set R will be equal to r-boundary nodes. We can 
then proceed to determine the degrees of n — f remaining nodes, or equivalently, partition the 
positive integer s (8) into n — f integers, each assuming a value between 1 and n — 1, and a 
lower bound for the degrees of r — f nodes [14]. The possible values for the partitioning comes 
from the fact that we expect the resulting graph be connected while respecting the bounds on 
the maximum allowable node degrees. 

Partitioning integers without constraints on the resulting partition is often referred to as 
unrestricted partitions. Restricted partitions, on the other hand, are those with constraints on 
the largest value of the partition that is no greater than a value of Ku, or no smaller than Kl, or 
both. Algorithms that generate unrestricted partitions can often be used to generate the restricted 
ones by certain modifications. Several such algorithms, dealing with unrestricted and restricted 
integer partitioning, have been suggested in literature. In the context of the graph realization 
using the proposed system identification method, we proceed to use the algorithms in [15], in 
order to generate different possible sets of n — r integers between 1 and n — 1, and r — f nodes 
with specified lower bounds on their degree such that their sum is s (8). 

C. Integer Partitioning Algorithms and Complexity Analysis 

Consider a degree sequence {di,d2, ■ ■ ■ , dn-f} while di + d2-\ — ■ + dn-r = s, with specified 
lower bound on r — f of them. Without loss of generality, assume that the first {di, . . . , dr-f} 
degrees are lower bounded as 




(8) 



di > Li for i — 1, . . . ,r — f. 



(9) 
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We are interested to find all possible partitioning of s into n — f integers between 1 and n — 1 
satisfying (9). We have the following observation; see [16]. 

Lemma III.1. Let the number of partitioning of s into n — f integers between 1 and n — 1 be 
denoted by Pn-f{s). Then 



The following algorithm, proposed in [15] finds all partitioning of s into m = n — f integers 
between 1 and n — 1 satisfying (9). The partitioning of s with m components can be generated 
in increasing lexicographic order by starting with di — d2 — ■ ■ ■ — d^-i = 1, dm — s — m + l 
and continuing as follows. To obtain the next partition from the first one, scan the elements 
from right to left, stopping at the right most di such that dm — di > 2. Replace dj by + 1 
for j = i,i + 1, . . . ,m — 1 and then replace dm hy s — Yl^=i '^j- example, if we have 
s = 12, m = 5, and the partition {1, 1, 3, 3, 4}, we find that 4 is greater by 2 than the rightmost 
1, and so the next partition is {1, 2, 2, 2, 5}. When no element of the partition differs from the 
last by more than 1, we are done. 

In the suggested algorithm the output size of each partitioning of s into some arbitrary number 
of integers m, Pm{s), is 0{s). This means that the total output size is 0{sP{s)). The approximate 
size of the number P{s) is provided by the following asymptotic formula. 



In other words, P{s) grows faster than any polynomial, but slower than any exponential 
function Q{s) = d^. However, in our application, we are interested in a subset of P{s) which has 
a specified size and satisfy certain constraints. Specifically, the integer s in (8) is approximately 
the number of edges in the graph. For simple graphs if 0{\£\) — 0{n), then the upper bound for 
the proposed partitioning is 0{ne^), and if 0{\£\) — 0{n^), the upper bound for the proposed 
partitioning is 0(n^e"). 

Although the algorithm above leads to a possible degree sequence for the underlying graph — 
consistent with the identification procedure- we need an additional set of conditions for ensuring 
that the obtained sequence in fact corresponds to that of a graph. 



Pn-f{.s) = Pn-f-l{s - 1) + (n - f)Pn-f{s - 1). 



(10) 
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Algorithm 1 Integer Partitioning 

di = (^2 = ■ ■ ■ = dm-i = 1, s-m + 1 

i = 1 

while i 7^ do 

if {di, d2, ■ ■ ■ , df} > {Li, L2, . . . , Lf} and I < di < n — 1, Vi then 

I output {di, da, - ^c^m} 
end 

i = m — 1 

while dm — di <2 do 

I i=i-l 

end 

if i 7^ then 

for j = m — 1 to i by -1 do 

dj = di + 1 

end 

end 

dm — S ^^j=i dj 

end 



Definition III.2. A graphical sequence is a list of nonnegative numbers that is the degree 
sequence of some simple graph. A simple graph with degree sequence d is said to realize 
d. 

Our next step is therefore to characterize the necessary and sufficient conditions for a set of 
integers to be graphical. For this, we resort to the following result. 

Theorem III.3. [17] For n > 1, an integer list d of size n is graphical if and only if d' is 
graphical, where d' is obtained from d by deleting its largest element A and subtracting Ifrom 
its A-th next largest elements. The only 1-element graphical sequence is di = {0}. 

Example III.4. Consider a sequence {3, 2, 2, 2, 2}. Since the number of odd degree nodes is 
odd, the sequence is not graphical. Let us also construct the d' sequences described above, as 
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{1, 1, 1, 2} and {1}, which again verifies that this sequence is not graphical. 

Since the maximum number of steps to check whether a sequence is graphical or not is n, 
the complexity of this algorithm is 0{ri). Given the particular algorithmic means of generating 
a graphical sequence for the required integer partitions, as detailed above, we now consider the 
problem of constructing graphs based on a graphical sequence. 

D. Degree Based Graph Construction Algorithms and Complexity Analysis 

Before describing the algorithm, we need to provide a few definitions.^ Let A{i) denote the 
adjacency set of node i defined as 

A{i) = {ttk I Ofe e V, afc > i, for all /c, 1 <k < di}. 

The reduced degree sequence d'|^(j) is obtained after removing node i with all its edges from 
Q. We now define the ordering < between two adjacency sets of node i, A{i) = {. . . , a^, • • •} 
and B{i) = {. . . , 6^, . . .}, as B{i) < A{i) if we have < for all 1 < /c < d^. In this case 
we also say that B{i) is "to the left" of A{i). The next lemma introduces a sufficient condition 
for the sequence d'|B(i) to be graphical. 

Lemma IIL5. [18] Let d = {di, d2, ■ ■ ■ , dn} be a non-increasing graphical sequence, and let 
A{i), B{i) be two adjacency sets for some node i e V, such that B{i) < A{i). If the degree 
sequence reduced by A{i) (that is d'|^(j)) is graphical, then the degree sequence reduced by 
B{i) (that is d'|B(j)) is also graphical. 

The above lemma guarantees preservation of "graphicality" for all adjacency sets to the left 
of a graphical one. Now consider a graphical degree sequence d on n nodes obtained from 
previously discussed integer partitioning approach. From the identity CAB — CAB, as we 
discussed in § 11, rxro/'^ entries of the system matrix A{Q) are known;^ define these set of edges 
as being "pre-determined" in the graph which cannot be repeated again. Put these connections 

^The necessary definitions and algoritiims have been discussed in [18] and are briefly described here to complement the 
presentation. 

*ln the case where C = , entries of CAB are known. 
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Algorithm 2 Degree Based Graph Construction 
Given a graphical sequence {di < ^2 < ■ ■ ■ < rf„ < 1}. 

I. Define the rightmost adjacency set ^^(i) containing the di largest index nodes different 
from i. Let us also define X{i) = {j eV,j^i s.t. {i,j} ^ £} as the forbidden neighbors 
of node i. Note that X(i) originally might contain some nodes according the forbidden set 
X{d). Create the set Ar(1) and X(l) for node 1 : connect node 1 to n (this never breaks 
graphicality). Set X{1) — {n}. Define the new sequence d' — {di — \X{l)\,d2, . . . , dn}. 
Let k — n — 1. 

LI. Connect another edge of 1 to k. Run the graphicality test in Theorem ELS. 

1.2. If this test fails, set = — 1; repeat I.l. 

1.3. If the test passes, keep (save) the connection, add the node k to the forbidden set 

and update the degree sequence d' = {di — \X{l)\,d2, . . . , dn}, set k = k — 1, and if 
i has edges left, repeat from I.l. 

II. Create the set ^(d) of all adjacency sets of node 1 that are colexicographically smaller 
than ^r(I) and preserve graphicality, i.e., 

^(d) = ^A{1) = {ai, . . . , ddj, e V I A{1) <cl Ar{1), d'U(i) is graphical} . 

where the operation <cl means a colexicographic order between two sets. 

III. For every A{1) G A{d) create all graphs from the corresponding graph reaUzation of d'|^(i) 
using this algorithm, where d'j^m is the sequence reduced by ^(1). 



in the forbidden set -'^(d). The Algorithm 2 describes how we can construct all possible graphs 
avoiding the edges in ^(d). When constructing A{d), checking graphicality is only needed for 
those adjacency sets which are incomparable by the ordering relationship to any of the current 
elements of A{d); for the remaining sets graphicality is guaranteed by Lemma IIL5. 

The total number of graphs that the above algorithm produces is Ili{di\). Of course, this 
procedure makes sense for degree sequences d for which there is only a small number of 
labeled graphs reaUzing it. An upper bound on the worst case complexity of the algorithm 
for constructing a sample from a given degree sequence d is Cd < 0{n\£\), with \£\ being the 
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number of edges in the graph. For simple connected graphs, the maximum possible number of 
edges is 0{n'^), and the minimum possible number is 0{n). If = 0{n), then Cd < 0{n'^), 

and if = 0{n'^), then Cd < 0{n^), which is an upper bound, independent of the degree 

sequence [19]. An algorithm to realize graphical degree sequences of directed graphs has been 
studied in [20] which can be use to extend these results to simple directed graphs. 

Our sieve method for confining our search for the underlying graph following the system 
identification of § II, thus involves: 

1) perform an integer partition on the value s (8), keeping the partitions that lead to a graphical 
sequence; let G denote the set of all graphs that remain after this first stage of the sieve. 

2) construct a candidate connected graph in G that is consistent with the matrix CAB and 
satisfy the given degree sequence, 

3) compare the Laplacian eigenvalues of the constructed graphs with the roots of the char- 
acteristic polynomial (7) and discard the inconsistent graphs. 

The candidate graphs for the underlying network topology are now among the ones that remain 
after this three step sieve. We note that the sieve method is guaranteed to reduce the search space 
for the graph structure by at least a factor of 2", a bound that is obtained from the bound on 
the number of permissible degree-based integer partitioning facilitated by the network system 
identification. An example for this procedure is given next. 

IV. An Example 

Our goal in this example is to gather information on the graph Q shown in Figure 2(a) using 
the system identification procedure. Using nodes 1, 2, and 3 as the input-output nodes in (2), 
we obtain (j)g{s) = + 220s^ + 190s^ + 804^3 + 1664s2 + 1344s. Since the polynomial 0g(s) 
has just one zero root, the underlying graph is connected. Moreover, the graph is not a tree due 
to the fact that a„_i 7^ 6. The graph has 11 edges and 224 spanning trees. We also obtain the 
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estimation matrices A, B, and C to be 



A = 



B = 



C = 



And also 



-2.0913 


9.5928 


-3.2434 


-1.9321 


-8.3119 


-3.6056 


-0.7817 


-2.6237 


-2.3217 


-3.3284 


-2.1037 


0.8747 


0.2080 


-0.1114 


-3.4715 


0.8443 


0.0218 


-0.2108 


-0.3304 


-0.1269 


0.5219 


-3.1500 


0.0999 


0.3707 


-0.6406 


2.3385 


-3.4079 


-4.7957 


-7.1871 


0.6395 


-1.8704 


9.1263 


-2.8730 


-1.3239 


-7.4520 


-3.4763 



-9.7903 13.7604 

1.6650 0.0982 

0.6348 0.1419 

0.4699 -0.3308 

0.3430 2.7623 



-4.0211 
-1.6752 
-0.3147 
-0.4116 

-3.1201 



-9.3019 12.8359 -3.5526 

-28.1542 2.3826 0.3547 
-28.2311 2.0005 0.0729 

-28.2455 2.5397 -0.5074 



0.8880 0.6434 30.0441 
0.2781 0.9634 30.1259 
-0.5056 0.5157 30.1419 



CAB 



-3.0000 
1.0000 
0.0000 



1.0000 -0.0000 
-4.0000 1.0000 
1.0000 -3.0000 



Since the diagonal of the matrix CAB is [—3, —4, —3]^, di = 3, d2 = 4, and = 3. In the 
meantime, the sum of degrees of the remaining nodes is 12. The possible integer partitions for the 
remaining three nodes such that the sum is 12 and each degree is less than 6 will be {5, 5, 2}, 
{5,4,3}, and {4,4,4}. Therefore, the set of the possible degrees sequences is comprised of 
{3, 4, 3, 5, 5, 2}, {3, 4, 3, 5, 4, 3}, and {3, 4, 3, 4, 4, 4}. According to the Theorem in.3, three set 
of integer partitioning are realization of some graphs. 

Next, we construct the graphs with the three candidate degree sequences by implementing 
the algorithm in [18] and constructing all connected graphs consistent with CAB , and with 
degree sequences {3, 4, 3, 5, 5, 2}, {3, 4, 3, 5, 4, 3}, and {3, 4, 3, 4, 4, 4}. The first degree sequence. 
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{3,4,3,5,5,2}, and the degree sequence {3,4,3,5,4,3} are realizations of two graphs while 
the third degree sequence is a realization of two graphs. All these graphs satisfy the constraint 
imposed by CAB. By comparing the Laplacian spectra for the corresponding 4 graphs with 
the roots of the identified characteristic polynomial, we can thereby identify the original graph. 
Three candidates for the constructed graph are depicted in Fig. 2(b)-(d). 

V. Conclusion 

In this paper, we introduced a network identification scheme which involves the excitation and 
observation of nodes running consensus-type coordination protocols. Starting with the number 
of vertices in the network as a known parameter, as well as the controllability and observability 
of the resulting steered- and-observed network, the proposed procedure strives to collect perti- 
nent information on the topology of the underlying graph. In this direction, we examined the 
applications of spectral characterization of graphs as well as a sieve method that is based on 
integer partitioning algorithms and feasible graphical sequences. 
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Fig. 1. The percentage of random planar graphs that are controllable from at least one node 
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(a) (b) (c) (d) 

Fig. 2. (a) a simple graph on 6 nodes, (b) a candidate graph constructed with degree sequence {3, 4, 3, 4, 4, 4}, (c) a candidate 
graph constructed with degree sequence {3, 4, 3, 5, 4, 3}, (d) a candidate graph constructed with degree sequence {3, 4, 3, 5, 5, 2}. 
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